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ABSTRACT 

Following the Chandra Orion Ultradeep Project (COUP) observation, we have studied the chemical 
composition of the hot plasma in a sample of 146 X-ray bright pre-main sequence stars in the Orion 
Nebula Cluster (ONC). We report measurements of individual element abundances for a subsample 
of 86 slightly-absorbed and bright X-ray sources, using low resolution X-ray spectra obtained from 
the Chandra ACIS instrument. The X-ray emission originates from a plasma with temperatures and 
elemental abundances very similar to those of active coronae in older stars. A clear pattern of abun- 
dances vs. First Ionization Potential (FIP) is evident if solar photosphcric abundances are assumed 
as reference. The results are validated by extensive simulations. The observed abundance distribu- 
tions are compatible with a single pattern of abundances for all stars, although a weak dependence 
on flare loop size may be present. The abundance of calcium is the only one which appears to vary 
substantially between stars, but this quantity is affected by relatively large uncertainties. 

The ensemble properties of the X-ray bright COUP sources confirm that the iron in the emitting 
plasma is underabundant with respect to both the solar composition and to the average stellar photo- 
spheric values. Comparison of the present plasma abundances with those of the stellar photospheres 
and those of the gaseous component of the nebula, indicates a good agreement for all the other ele- 
ments with available measurements, and in particular for the high-FIP elements (Ne, Ar, O, and S) 
and for the low- FIP element Si. We conclude that there is evidence of a significant chemical fraction- 
ation effect only for iron, which appears to be depleted by a factor 1.5-3 with respect to the stellar 
composition. 

Subject headings: stars: cluster — stars: activity — stars: coronae — stars: late-type — X-rays: stars 



1. INTRODUCTION 

Chemical composition is one of the key properties of 
astrophysical environments, fundamental for classifying 
stellar populations and studying the evolutionary history 
of galactic chemistry over different spatial scales. Sev- 
eral processes altering clement abundances in the vicinity 
of individual stars are particularly efficient in the early 
phases of stellar evolution: selective trapping in grains, 
high-energy photon and particle irradiation of the cir- 
cumstellar medium, mass exchange between stars and 
protoplanetary disks via accretion and outflows, and frac- 
tionation effects in stellar coronae and magnetospheres. 
In pre-main sequence systems, these physical processes 
occur on short time scales and cause dramatic changes 
as new planetary systems form out of the circumstellar 
nebula. One of the widely discussed issues is why planet- 
hosting stars appear to be characterized by a mctallicity 
hig her (by 0.24 dex, on average) than stars without plan- 
ets ([Santos et al.ll2005l and references therein). 

The Orion Nebula Cluster (ONC) is one of the best 
studied star forming regions in the sky, and the chem- 
ical composition of the associated H II region has been 
historically considered a st andard reference for ionized 
gas in the nearby Galaxy (Es teban et al.l 120041 ) . Op- 
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tical spectroscopy of the nebula is the traditional ap- 
proach employed to study the abundances of important 
elements in this environment, as well as in other star 
forming regions, because of the difficulty in obtaining 
photospheric abundance measurements for faint, rapidly 
rotating young stars. Abundance studies of Orion stars 
have been performed ma i nly o n B m a in-sequence mem- 
bers (ICunha fe Lambert! fl99l 1 1994 ISimon-Diaz et all 
120061 : ICunha. Hubenv. fe Lanzll2006l ). Only a few mea- 
surem ents are available for slowly-ro tating F and G 
stars (ICunha. Smith, fc LambertJll998f) and K-M mem- 
bers (|Cunha fc Smithll2005h . ~ 

An alternative way to determine the chemical com- 
position o f late-type s tars is emerging from X-ray spec- 
troscopy (|Giidell (200l Favata & Micela 2003, and ref- 
erences therein). The ONC was selected in 2003 for a 
very long observation program in X-rays with the Chan- 
dra satellite, kn own as the Chandr a Orion Ultradeep 
Project (COUP. iGetman et al.l [20051 ). This program is 
providing an unprecedented wealth of information about 
the stellar population of the ONC and various character- 
istics of this prototypical stellar and planetary nursery. 

Among the s alient global properties of the ONC, 
iFeigelson etaLl (|2005D noted the presence of a strong 
spectral feature around 1 keV in the cumulative spec- 
trum of all detected X-ray sources, identified with the 
emission line complex due to H-like and He-like Ne ions 
in hot plasma associated with ~ 1400 ONC members. 
The prominence of this feature suggests a high abun- 
dance of Ne in the X-ray emitting plasma, a character- 
istic already observed in other magnetically active stars 
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(see review by iGiidell l2004f) together with an apparent 
depletion of iron in corona with respect to the expected 
photosphcric composition. 

This behavior is linked to the First Ionization Potential 
(FIP) of the elements and possibly to the stellar activ- 
ity level: the composition of the coronal plasma in the 
Sun, and particularly in long-lived magnetic structures, 
appears enriched by low-FIP elements (FIP < 10 eV) 
by about a factor 4 with respect to photospheric values. 
This is known as the FIP effect (e.g. iFeldman fc Laming) 
2000), a characteristic observe d also in other low-activity 
stars (jFavata fc Micelal 120031 . and references therein). 
High-activity RS CVn-type and Algol-type binaries stars 
exhibit a different behavior with a tendency for low-FIP 
elements such as iron to become depleted with respect 
to high-FIP elements like argon and neon. This is called 
"inverse FIP effect" (|Brinkman et al.ll200lD . 

The above scenario requires further investigations for 
several reasons. First, photospheric abundances are usu- 
ally uncertain due to severe NLTE and rotational broad- 
ening effects on optical spectra of active stars, and solar 
values are often employed as the reference for the stellar 
coronal abundances. The solar photosp heric compositio n 
itself has been recently questioned by lAsplundl (|2005f) . 
based on detailed 3-D modeling of the solar atmosphere. 
Moreover, the photospheric Ne abundance is not directly 
measurable even in the Sun because no Ne lines occur at 
optical wavelengths. 

Second, the driving mechanism(s) for FlP-related frac- 
tionation in stellar upper atmospheres still escapes a 
clear understand i ng alt h ough some models have emerged 
(lArge fc Mullanl 119981 : ISchwadron. Fisk. fc Zurbuchenl 
19991: lLaming||2004h . 

Third, the situation is made more complex by 
the presence of very high Ne abundance ratios in 
two classical T Tauri stars (CTTS), TW Hya and 
BP Tau, where the soft X-ray emission is often at- 
tributed to an accretion shock rather than to magnetic 
activity (IKastner et al.l 120021 : IStelzer fc Schmittl 12004 
iSchmitt et al.l I2005D . An alternative explanation for 
these high Ne abundances has been suggested based on 
an origin of the plasma in gas accreted from the cir- 
cumstellar disks where refractory elements may be de- 
plet ed into solids undergoing growth into pl anetary bod- 
ies (|Drake. Testa. fcHartmannf2005)l l2005h . However, 
X-ray variability characteristics strongly favors CTTS X- 
ray emissi on dominated by magnetic flares rathe r than 
accretion (|Stassun et al.ll2006t IStelzer et al.ll2007f l. 

The COUP observation provides us with the largest ho- 
mogeneous sample ever studied to address the issue of el- 
ement abundances in X-ray emitting plasmas associated 
to young stars. Here we can exploit ensemble statistical 
properties to overcome uncertainties in individual stellar 
measurements. The uncertainties on abundances derived 
from X-ray spectra, even those obtained with the high- 
est spectral resolut ion, may be larger th an formal sta- 
tistical error bars (Schmitt fc Nessl 120041 : iMaggio et alJ 
l2005f ) . The large COUP stellar sample can compensate 
for this difficulty. 

To address these issues, we present in this paper the 
results of a detailed analysis of CCD-resolution X-ray 
spectra of ONC members to derive the abundances of in- 
dividual elements for a large number of young stars. In 
fj2j we introduce the observation and the sample selec- 



tion, fj3]is devoted to the methodology of analysis, while 
the results are presented in 2] and discussed in §S}5}-9. 

2. OBSERVATION AND SAMPLE SELECTION 

The COUP data were obtained in January 2003 
with the Advanced CCD Imaging Spectrometer (ACIS; 
iGarmire et al.ll2003f) on-board the Chandra X-ray Obser- 
vatorv (jWeisskopf et al.|[2002f ) by combining six consec- 
utive observations of the Orion Nebula Cluster (ONC) 
with the same aimpoint and roll angle. The result- 
ing data set has a total exposure time of ~ 838 ks (9.7 
days) distributed over 13.2 days. The data reduction 
and analysis resulted in 1616 detected X-ray sources, 
of which about 140 8 are associated with ONC stars 
(jGetman et al.|[2~005l ) . X-ray spectra and light curves for 
each source were constructed from events collected in a 
polygonal region centered at the source position includ- 
ing ~ 90% of the encircled energy. Here we employ the 
spectra and the related auxiliary response files (ARFs), 
created with the ACIS_extract software, as explained 
in Getman et al. 

The criteria for sample selection were motivated by the 
need to perform a detailed analysis of the ACIS spectra, 
and - more specifically - to determine both the plasma 
temperature distribution in the corona and the abun- 
dances of a number of individual elements. The main 
requirement for this kind of analysis is a very strong sig- 
nal (noise is negligible in these cases) and an excellent 
knowledge of the instrument spectral response. We se- 
lect stars with at least 5000 total (net) counts in the 
0.5-8 keV energy band. Source s heavily affected by pile- 
up (flagged "a" in Table 6 of iGetman et alJ 120051 ) are 
excluded but the sample includes 23 stars with "mild" 
pile-up effects ("w" flag) for which the usual whole ex- 
traction region was used. We exclude four early-type 
stars (T c ff > 10, 000 K) as their X-ray emission likely 
originates in wind shocks rather than in a corona. How- 
ever, we have retained 34 stars with unknown T e gp. The 
final sample contains 146 COUP sources listed in Table 
[TJ about half of them have more than 8000 total counts, 
and 55 exceed 10 4 counts. 

Further a posteriori selection was made based on the 
amount of interstellar absorption exhibited in the X- 
ray spectra. Reliable determination of the abundances 
of some elements is possible only for slightly absorbed 
sources. We will focus our attention on the 86 sources 
with intervening hydrogen column densities Nh < 6 x 
10 21 cm~ 2 . 

The optical properties of the 146 stars in our initial 
sample are reported in Table [TJ based on the tables in 
IGetman et all (|2005l ). Fi gure [T] shows two near-infrared 
(NIR) color-color diagrams for all the COUP detected 
sources with available photometric measurements in the 
NIR bands of interest. These diagrams show the theo- 
retical locus of stars at the appropriate ONC age, and 
loci populated by objects with different amounts of red- 
dening. Twelve sources fall to the right and below these 
loci in both color-color diagrams indicating they prob- 
ably have dusty circumstellar disks in addition to red- 
dening. Three of these sources plus two more without 
NIR photometric excesses (COUP 382, 579, 597, 758, 
and 1409) reveal protoplanetary disks seen in silhouette 
against the bright nebula ( "proplyds" ) in Hubble Space 
Telescope imaging (iKastner et al.l 120051) . 
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Fig. 1. — Near-infrared color-color diagrams of the COUP 
sources with the stars in the present study indicated by filled red 
circles. The number of stars in the two diagrams differ because of 
incomplete IR photometric data f or some of them. In magenta, the 
locus of stars with ages 1-5 Myr ( Sics s. Dufour. fe Foresti"n"ill2000l , 
using T e g-color conversions by Kcnyon & Hartmann 1995). The 
black arrow indicates a reddening vector corresponding to A v = 10, 
the dashed green lines bracket the region where stars with no NIR 
excess are expected, and the b lue segment marks the locus of dered- 
dened classical T Tauri stars (Meyer et al. 1997). 

Table [T] also reports the eq uivalent width of th e 8542 
A Ca II line from the study of iHillenbrandl (|1997t) . In 14 
cases, this line is in emission with equivalent width larger 
than 1 A, suggesting that the central star is actively ac- 
creting material from the circumstellar medium. 

3. ANALYSIS 

Essentially all X-ray bright COUP sources show 
significant var i ability of their X- r ay emission level 
(IGetman et al.l 120051 iFavata et ail l2005t IWolk et all 
I2005F In~ many cases, this variability is associated 
with large flares which are well-fit by a solar-type flare 
model where plasma is suddenly heated by a magnetic 
reconnec tion event and co ols on timescales of hours- 
to-days (|Favata et alJ [2005) . In the present work, we 
perform the X-ray spectral analysis collecting photons 
over the whole observation length which includes qui- 
escent and flaring episodes. This may confuse inter- 
pretation of elemental abundances, as both tempera- 
tures, and sometimes FIP effects, vary during stellar 
flares (e.g. lAudard. Giidel. fe, Mewel 120011 : lOsten et all 
120041: iNordon. Behar. fe Giidel I2006T ) . Our results thus 
give thermal and chemical characteristics of the emitting 
plasma which are spatially and temporally averaged over 
for each star. This conflation of 'quiescent' and flare con- 
ditions may not be avoidable as there may not exist any 
true 'quiescent' state in very active stars such as those in 
the ONC because, even during periods of little variabil- 
ity, the X-ray emission likely aris es from a sup erposition 
of a multitude of weaker flares (jGudell 120041 and refer- 
ences therein) . This continuous flaring paradigm implies 
that the abundance properties we derive for the COUP 
sources probably describe time-averaged dynamical con- 
ditions rather than static equilibrium conditions of the 
X-ray emitting plasma. 

We assume that the observed emission can be mod- 
eled as a collisionally-excited plasma in ionization equi- 
librium, and we adopt the emissivities predicted by the 
Astrophysical Pla sma Emission Code (APEC VI. 3.1, 
iSmith et al.l l200lh in the spectral fitting p rocess. This 
choic e, rather than the MEKAL emissivities ijMewe et al.l 
119951 ) adopted in previous COUP works, is motivated by 
the significantly larger number of emission lines and more 
updated atomic data of APEC vs. MEKAL as imple- 
mented in the XSPEC spectral analysis package. This is 



especially important in our work, aimed to derive infor- 
mation from line complexes due to specific atomic species 
of individual elements, which can be resolved at most 
only marginally in the available CCD spectra. In any 
case, this choice may affect the abundance measurements 
of some elements, but not our global results. 

The data sets and instrument spectral response for 
each sourc e are obtained from the data reduction de- 
scribed bv IGetman et al.l (|2005l ). Special features of the 
COUP data processing include a 0.5 — 8 keV energy range 
and, for these strong sources, grouping of energy chan- 
nels such that there are at least 60 net counts in each 
spectral bin. 

3.1. Spectral diagnostics of element abundances 

We adopt a global spectral fitting approach with multi- 
temperature models where individual element abun- 
dances are free parameters, in addition to the tem- 
perature and volume emission measure of each compo- 
nent and the interstellar hydrogen column density Ah 
to the star. Photoelectric absorption in the interstellar 
m edium (ISM) is modeled with c ross-sections obtained 
bv lMorrison fe McCammonl (|1983h . The spectral resolu- 
tion of ACIS CCDs is such that a model with two or three 
isothermal components usually provides a good descrip- 
tion of the observed spectra. This modeling approach 
provides an approximate description of the continuous 
distribution of temperatures undoubtedly present in the 
X-ray emitting plasma, but it is certainly adequate to 
the amount of information provided by instruments with 
resolution power w 10-30. 

For the abundance measurements, we consider the el- 
ements O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni. At 
the coronal temperatures typical of young active stars 
(T w 10 7 K), all of these atoms have important H-likc 
or He-like ion lines in the ACIS wavelength range (1.5- 
27.6 A, Fig. [2]). Iron and Nickel are also represented by 
a large number of L-shell emission lines. However, clear 
spectral signatures of each element depends on several 
factors: relative line emissivities, plasma emission mea- 
sure vs. temperature distribution, line blending, inter- 
stellar absorption and, of course, the abundance of each 
element in the plasma. Due to the complexity of these 
dependencies, we have performed extensive sets of simu- 
lations to validate our spectral fitting results. Details on 
these simulations are reported in Appendix and will be 
included in discussion of specific results below. 

The neon abundance is especially interesting for sev- 
eral reasons. First, this abundance cannot be determined 
in stellar photospheres because of the lack of suitable 
absorption lines in optical spectra, hence X-ray data 
provide the best opportunity to perform such a mea- 
surement. Second, neon has the highest First Ioniza- 
tion Potential (FIP = 21.56 eV) of any atoms except he- 
lium, hence its abundance with respect to iron (with a 
low FIP) provides crucial information about the strati- 
fication of elements with different FIPs in stellar atmo- 
spheres (§HJ. Third, H-like and He-like Ne ions produce 
prominent emission lines at energies around 1 keV where 
Chandra/ ACIS sensitivity is highest. The correct deter- 
mination of the Ne abundance is not free of difficulties 
due to the proximity of its most intense emission lines 
with L-shell iron and nickel lines. Figure [3] shows this 
with a plot of the integrated emissivity of Ne, Mg, Fe, 
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Fig. 2. — Examples of spectra of Orion COUP sources with 
major emission line complexes labeled by their element abundances 
(in solar photospheric units). The best-fit 3-T model parameters 
(Nn in units of 10 22 cm -2 , T in keV, fc in cm -5 ) and the reduced 
X 2 are indicated in the top panels, while bottom panels show the 
residuals in units of standard deviation. 



and Ni lines in the wavelength range 9-14.5 A, vs. plasma 
temperature for a solar mixture of chemical elements 
(|Anders &; Grevesse] 1 1989ft . Although the Fe emissivity 
exceeds that of Ne for solar abundances, when the Ne/Fe 
ratio is sev eral times the sola r ratio as foun d in many ac- 
tive stars ([Drake et al.ll200H iGtideH 12004). the strength 
of the Ne lines in this spectral range becomes compara- 
ble to or greater than that of the iron lines. Thus, Ne 
abundances in active stellar coronae are sufficiently high 
to bring this element within the reach of global spectral 
analysis, as we will illustrate shortly. 
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Fig. 3. — Total emissivity of emission lines from Ne, Mg, Fc, 
and Ni ions in the wavelength range 9— 14.5A, vs. temperature from 
th e APEC plasma model as suming the solar element abundances 
of Endcrs & Grcvcsse] lfl989Tl . 



3.2. Spectral fitting procedure 

Our procedure starts with fitting of 2-temperature (2- 
T) and 3-temperature (3-T) plasma models to all the 
source spectra using a x -minimization a lgorithm imple - 
mented in the XSPEC (VI 1.3) package (|Arnaudlll996l) . 
This step is repeated a few times with different start- 
ing values of the free parameters to avoid local minimum 
X 2 solutions. An F-test is then applied to determine 
whether the (usually lower) x 2 obtained with the 3-T 
model represent a significantly better fit to the data, or 
rather the improvement is entirely due to the larger num- 
ber of free parameters with respect to the 2-T model. We 
adopt 3-T results only if the following two conditions are 
met: (i) the 2-T model has a poor fit with probability 
P{X 2 ) < 10%; and (ii) the F-test shows the 3-T model 
is significantly better with P(F) < 10%. The unneces- 
sary introduction of a third thermal component, which 
typically has the lowest temperature, may alter the abun- 
dances of elements (such as O, Ne, and Mg) with emission 
lines only in the low-energy tail of the spectrum. 

Our elemental abundances are s caled to the widely 
used solar system abundances of lAnders k, Grevessi 
(fl989h . In § El we will also discuss the imp lications of the 
recent solar composition recommended bv lAsplund et all 

This procedure resulted in 118 spectra fitted with a 
2-T model and 28 spectra fitted with 3-T models. For 
18 spectra, none of the models provides a statistically 
acceptable fit at the 99% confidence level, and in 12 
of these cases the 3-T model is not significantly better 
than the 2-T model. This suggests that the poor fit is 
not due to the limited number of components adopted, 
but rather to other causes, perhaps residual problems in 
the calibration of the instrument response. Inspection 
reveals in most of these cases large residuals near the 
iridium edges at 5.7-5.9A associated with the coating 
of the Chandra mirror. A formally better x 2 could be 
obtained by ignoring a narrow wavelength range in this 
spectral region, with little variations of the best-fit pa- 
rameters. For several sources, we cannot exclude also an 
astrophysical process such as a departure of the plasma 
conditions from thermal equilibrium or temporal and/or 
spatial variations of element abundances in the emitting 
plasma. Recall that most of the sources are strongly 
variable and characterized by important flaring events. 
Nonetheless, since source variability is the norm rather 
than exception for ONC stars, we have not discarded any 
source based on a variability criterion. 

3.3. Validation of spectral parameters 

We have estimated the uncertainty on individual quan- 
tities derived from these highly nonlinear spectral fits 
using the XSPEC code, with the criterion A% 2 = 2.7 
around x 2 minimum, corresponding to the 90% confi- 
dence level for one int eresting model parameter at a time 
(|Lampton et al.|[l976f ). In particular, we have evaluated 
uncertainties in Fe and Ne abundances, elements with 
the strongest line signatures in our spectra, by allow- 
ing temperatures, emission measures, and abundances of 
the four elements with lines in the wavelength range 9- 
14.5 A, (Ne, Mg, Fe, and Ni; see Figured]) to vary freely. 
The resulting uncertainties are mostly less than a factor 
of two around the fitted value; more precisely, for half 
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Fig. 4. — Left: Comparison of absorption column densities ob- 
tained by iGetman et al, (2005) with those in the present work. 
Right: Plot of absorption column density from COUP spectral fits 
vs. visual absorption. The two curves give the gas-to-dust relation- 
ship Nu = 1-6 X 10 21 A V (lower) and Nh = 2.2 X 10 21 A V (upper). 
Full symbols are used for stars in the count-limited sample. 



the measurements the relative error is within 30%, and 
it exceeds a factor 2 only in 5% of the cases (see error 
bars plotted in Figure [5]) . 

We have evaluated the reliability of the derived abun- 
dances and their uncertainties in the simulations de- 
scribed in Appendix A. We find that the overall pattern 
of abundances for most elements is recovered with little 
bias by our analysis procedure, although some elements 
(Ca, Ni, Mg and O in particular) could be vulnerable 
to systematic errors. The XSPEC errors for Fe and Ne 
abundances are sometimes smaller than indicated by the 
simulations (§ I4.2[) . 

One of the sources of uncertainty could be the ac- 
tual emission measure distribution vs. temperature in 
the emitting plasma. To test the possibility that our 
element abundances derived from global spectral fitting 
could be affected by inadequate X-ray emission models, 
we have performed simulations with input emission mea- 
sure distributions more complex than simple 2-T or 3-T 
approximations. We then checked that the plasma abun- 
dances are correctly recovered in spite of the mismatch 
between the actual source emission measure distribution 
and the adopted fitting model. For a few sample stars, 
we also fitted the observed X-ray spectra with alterna- 
tive plasma emission models, having a fixed grid of tem- 
peratures and variable emission measures, and we have 
obtained abundance measurements consistent with the 
results presented above, within statistical uncertainties. 
These simulations and tests give us confidence that our 
spectral analysis provides us with a reasonably accurate 
description of the abundance patterns in the observed 
coronal plasma, at least for the sources without strong 
interstellar absorption (see below). 

A final check on our spectral fitting is shown in Figure 2] 
where the derived interstellar hydrogen column densities 
(Nh) are compared with values obtained by Getman et 
al. (2005) for all COUP sources and with dust redden- 
ing estimated from optical spectroscopy. Our values are 
very closely correlated to the COUP values except for 
a systemic offset by about 0.1 dex. We attribute this 
small discrepancy to the differences in the spectral fitting 
procedure used in the COUP analysis (1-T/2-T MEKAL 
plasma models rather than 2-T/3-T APEC models). The 
scatter in the Ah — Ay plot is similar to that seen in the 
full COUP sample. 

As the Ah values in our sample exhibit a wide range 
from ~ 2 x 10 20 to ~ 10 23 cm~ 2 , we have investigated 
whether absorption might affect the elemental abun- 
dance measurements. Figure [5] shows a strong increase 
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Fig. 5. — Scatter plots of best-fit abundance measurements 
of Fe (left) and Ne (right) vs. H column density. Plus and arrow 
symbols indicate upper limits. The heavy vertical line is drawn at 
the threshold value we have adopted to define our low-absorption 
sample. 

in the spread of Ne and Fe abundances, accompanied by 
larger statistical errors, as Ah increases. To minimize 
the influence of the above effects but still keeping the 
sample size as large as possible, we apply a threshold 
A H < 6 x 10 21 cm~ 2 (corresponding to Ay < 3-4) to 
avoid the high scatter in Ne and other abundances due 
to absorption. At this threshold, the attenuation at the 
Ne X Lya wavelength (12.13 A) is about a factor 4. There 
are 86 sources in our low-absorption sample, and we will 
focus our attention on them in the next sections. 

4. RESULTS 

The best-fit spectral model parameters for the 86 stars 
in the low-absorption sample are reported in Table [2] A 
subset of 35 sources in this sample have more than 10 4 
counts in their spectra, and we will call it the count- 
limited subsample. The table gives the derived absorp- 
tion, plasma temperatures and emission measures, abun- 
dances for 9 elements, the reduced x 2 °f t ne fit, and 
the source X-ray flux in the 2-8 keV band. The me- 
dian values of the abundance distributions and the cen- 
tral 68% ranges are reported in Table [3l for both the 
low-absorption sample and the count-limited subsample. 

4.1. Coronal temperatures and elemental abundances 

Figure [6] shows boxplots of temperatures, ratios of 
emission measures, and H column densities for the low- 
absorption subsample. These ONC stars are character- 
ized by coronal plasma with temperatures ranging from 
w 5MK to 25 MK (median values for the 2-T or 3-T 
models). The high-temperature components are domi- 
nant in most cases with emission measures typically two 
times larger than for the cool (T < 10 MK) components. 
It is worth noting that the thermal characteristics of our 
sample stars are optimal for measuring Ne abundances, 
together with those of Mg and Ni, because the emis- 
sivities of the relevant emission lines peak in the same 
temperature range (see Figure [3|). 
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Fig. 6. — Box plots of temperatures, ratios of emission measures, 
and H column densities, derived from 2-T and 3-T fits (see Table 
[2j. The upper and lower edges of each box comprise the central 
68% of the data, the central value is the median. Squares mark 
individual measurements. 
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Fig. 7. — Box plots of the best-fit abundance distributions for 
each element, sorted by increasing First Ionization Potential. The 
boxplot on the left in each pair refers to the full sample of slightly 
absorbed sources 

(JV H < 6 X 10 21 cm -2 ), while the boxplot on the 
right pertains to the count-limited subsample (> 10, 000 extracted 
counts). 

Boxplots of best-fit abundance values vs. First Ioniza- 
tion Potential (FIP) for the low-absorption sample and 
for the count-limited subsample are shown in Fig. [7l 
The three values indicated by each box (lower and upper 
edges, and central segment) represent the 68% range and 
the median reported in Table [3J 

The striking feature of all these plots is the systematic 
pattern of abundance values vs. FIP: relatively low abun- 
dances with respect to the solar photosphcric composi- 
tion are consistently found for low-FIP elements Fe, Mg, 
and Si while elements with higher FIP are increasingly 
more abundant. Calcium and nickel abundances do not 
follow this trend, and our simulations (Appendix A) con- 
firm the reliability of this result, in spite of some possible 
systematic error in the case of the Ca, and the sensitivity 
to line blending effects of the Ni measurement. 

4.2. Reliability of the spectral analysis results 

Before attempting any interpretation of the results we 
need to discuss their robustness against a number of 
possible sources of uncertainty. We first considered the 
count-limited subsample comprising the 35 sources with 
more than 10,000 counts, and the subsamples of sources 
with 2-T or 3-T best-fit models (74 and 12 sources, 
respectively). These yield essentially the same abun- 
dance distributions as the low-absorption sample, al- 
though with slightly different amount of scatter. Hence, 
the results do not depend on the source strength or as- 
sumed plasma temperature distribution. 

We performed several simulations as described in Ap- 
pendix A to investigate a variety of other possible ef- 
fects. The simulations provide us with distributions of 
the best-fit abundances which take into account: the pho- 
ton counting statistics at each wavelength, the possible 
cross-talk between elements with emission lines falling 
within the instrument spectral resolution, and the cross- 
talk between line strength and continuum level deter- 
mined by the normalization of the thermal components. 
These simulations show that our procedures reliably re- 
cover true source coronal abundances. The observed FIP 
pattern is recognized self-consistently only in simulations 
with input models assuming that specific pattern, and 
FIP effects are not artificially introduced when they are 
not present. 

Uncertainties of the model abundances are evaluated 
through the simulations and found to have scatter sim- 
ilar to that seen in the observed sample. For 2-T or 
3-T source model spectra and count rates typical of our 



COUP sources, our spectral analysis is able to recover 
the input values within a factor of 2 for Fe, Si, and S, 
and within a factor of 3 for Ni, O, Ar, and Ne. Mg and 
Ca abundances are the most uncertain by factors 5-10. 

5. CORONAL ABUNDANCES IN THE X-RAY 
LUMINOUS ONC STARS 

It is important to recognize that the sample of 86 ONC 
stars giving the results in Tableland Figure[7]is uniquely 
large and homogeneous in the field of stellar X-ray spec- 
troscopy. Abundance measurements based on ACIS low- 
resolution CCD spectra are usually impractical due to 
either insufficient counts or pileup effects in high count 
rate stars. Abundance measurements based on high- 
resolution grating X-ray spectroscopy with Chandra and 
XMM-Newton are available up to now for less than 30 
late-type stars with vastly different ages and in disparate 
astrophysical environments. 

Given our large sample size, we can focus our atten- 
tion on the subsample of 35 sources with more than 
10,000 extracted counts which provides the most reli- 
able abundance measurements. We have verified using 
Kolmogorov-Smirnov tests that the elemental abundance 
distributions for this subsample, are statistically indistin- 
guishable from the distributions obtained for the sources 
having between 5000 and 10000 counts. However, our 
entire study is certainly biased toward the most magneti- 
cally active X-ray stars and is limited to slightly absorbed 
COUP sources. The characteristics of the more embed- 
ded Orion stellar population and of the stellar coronae 
with relatively lower X-ray luminosities are not treated 
here. 

The next step is to establish whether our sample of 35 
bright sources is consistent with a single distribution of 
coronal abundances without star-to-star variations. Fig- 
ure [8] shows a comparison between the observed and sim- 
ulated distributions of abundances for Fe and Ne. The 
simulation here was performed assuming a 2-T model 
with all parameters fixed at the median values of the 
observed distributions and taking into account the ac- 
tual photon counting statistics of the observed spectra 
(Appendix A). A Kolmogorov-Smirnov test finds that 
the observed distributions are consistent with the sim- 
ulations (P ~ 10 — 20%). Similar results are obtained 
for all the other elements, except for the peculiar case of 
the Ca (see below). Thus, the observed spread of abun- 
dances for each element is compatible with being due 
to the uncertainties on the measurements. There is no 
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Fig. 8. — Cumulative distributions of Fe (left panel) and Ne 
(right panel) abundances: results from simulations employing 2-T 
models (dotted line) are compared with 2-T spectral fitting results 
for our count-limited subsample (heavy dashed line) , and with the 
results for the full sample of slightly-absorbed sources (solid line). 
For both elements, the observed and simulated distributions are 
statistically undistinguishable at the 99% confidence level. 
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Fig. 9. — Left: Cumulative distributions of the Ne/Fe abundance 
ratios for the full sample of slightly-absorbed sources (solid line) 
and for the count-limited subsample (heavy dashed line). Right: 
Scatter plot of Ne vs. Fe abun danc es, with error bars evaluated 
from individual spectral fits (§ 13.20 . for the stars in the count- 
limited sample. The two straight lines indicate different Ne/Fe 
abundance ratios. 

clear evidence for different coronal compositions among 
the ONC stars in the count-limited sample. However, 
we show below that marginally significant differences be- 
tween certain subsamples may be present. 

Inspection of Table [3] shows that the iron abundance 
of the coronal plasma in our full sample of X-ray bright 
ONC stars is well constrained in the range 0.12-0.37 
times the solar value. The abundances of the low-FIP 
elements Mg and Si are compatible with the iron abun- 
dance, while the higher FIP elements S, O, Ar, and Ne 
appear systematically higher. Considering Fe and Ne as 
representative of the low-FIP and high-FIP species, re- 
spectively, we find a median Ne/Fe abundance ratio of 
6±2 (Fig. [9] left). While the uncertainties on the mea- 
surements for individual stars can be large (Fig. [9] right), 
the median value is well-measured and the possibility 
that Ne/Fe=l is confidently excluded. 

Nickel and calcium abundances do not follow a simple 
FlP-abundance relationship. The relatively high abun- 
dance of the low-FIP Ni may appear suspicious. The 
atomic database we have employed contains a large num- 
ber (> 1600) of L-shell Ni lines in the range 5-24 A, 
produced by all ions from Ne-like Ni XIX to Li-like 
Ni XXVI, which form at temperatures ranging from 8 MK 
to 25 MK, hence it is quite complete in this respect. The 
most prominent spectral lines are those of Ni XIX at 
12.44, 12.66, 13.78, 14.04, and 14.08 A, which fall close 
to the important H-like and He- like Ne lines, and a cross- 
talk between the two abundance parameters is possible. 
However, our simulations indicate that a high Ni abun- 
dance can be correctly recovered (Appendix A). We con- 
clude that the high Ni abundance looks real, although 
the uncertainties may be larger than for other elements. 




Fe abundance Ne abundance 

Fig. 10. — Cumulative distributions of Fe (left panel) and Ne 
(right panel) abundances for the count-limited sample (solid line) 
and for three different stellar groups: X-ray sources associated with 
proplyds (dotted line), stars with Ca II in emission (heavy dashed 
line), and stars with evidence of large flares (dot-dashed line). 

The best-fit calcium abundances are zero for about 



70% of the stars in our sample, independently from the 
amount of hot plasma indicated by the best-fit model. 
But Ca K-shell lines are clearly visible at ~ 3 A in the 
spectrum of several COUP sources and give high mea- 
sured Ca abundance values (Table [2|). Other important 
L-shell lines from Ca XVI- XVIII fall in the range 19- 
24A, i.e. at the end of the inspected ACIS band; these 
lines are much weaker than the O VII- VIII lines occur- 
ring in the same spectral region, and hence they are not 
useful for the determination of the Ca abundance. Sim- 
ulations performed with the Ca abundance set to zero 
predict a distribution of best-fit Ca values which is be- 
low the observed distribution, and simulations assuming 
a solar Ca abundance give a higher distribution (Ap- 
pendix A) . This suggests either an intrinsic spread in Ca 
abundance is present in the sample, or that some un- 
known systematic error in the analysis affects calcium 
abundance estimates. At present, we believe that the re- 
ported underabundance of Ca in most sources is a solid 
result. 

6. ABUNDANCES AS A FUNCTION OF OTHER 
PROPERTIES 

The above analysis indicates that the stars in our sam- 
ple, chosen to have very high X-ray luminosities (L x in 
the range 10 29 8 -10 31 erg s _1 ), share similar temperature 
distributions and chemical abundances. These similari- 
ties represent a major result of the work presented here 
and suggest that a single physical mechanism is operative 
in the sample. 

We can nonetheless investigate whether interesting 
subsamples behave as the whole population of X-ray- 
bright ONC stars. Inspection of the optical character- 
istics of our ONC sample (Table QJ reveals that our sam- 
ple includes nine stars with Ca II in emission suggesting 
active accretion (COUP 11, 66, 112, 141, 567, 579, 670, 
801, and 1608), and five stars associated with imaged 
proplyds (COUP 382, 579, 597, 758, and 1409). Sev- 
enteen stars in our sample were studied by Favata et al. 
(2005) for the presence of large flares in their COUP light 
curves; six (COUP 43, 141, 669, 752, 848, 1608) showed 
evidence of very long (L > 5i?*) flaring magnetic loops. 

In Figure[l0l we compare the Fe and Ne abundance dis- 
tributions for our count-limited sample with abundances 
of these three groups. We find that the X-ray sources 
associated with proplyds show on average higher abun- 
dances with respect to the sources in the count-limited 
sample, while the strong-Ca II stars and the stars in 
the flaring group are indistinguishable from the full sam- 
ple. However, even in the former case, the distributions 
for the stars in the subsample are not significantly dif- 
ferent from those in the count-limited sample at 90% 
Kolmogorov-Smirnov confidence levels. 

The case of the stars caught during strong flares is 
particularly interesting because time-resolved analyses of 
large flaring events in active stars have indicated in many 
cases an apparent increase of the plasma metallicity at 
the onset of the flare w ith respect to the quiescent phase 
(|Favata fe Micelall2003l . and references therein). We thus 
might expect that the COUP stars whose X-ray emission 
is dominated by strong flares could have higher Fe and Ne 
abundances. We do not observe this effect in Figure flOl 
where the flaring stars are considered all together. 

A more intriguing case is offered by Fig. 1 1 X I which 
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Fig. 11. — Cumulative distributions of Fe (left panel) and Ne 
(right panel) abundances for the count-limited samp le (solid line) 
and fo r two subgroups of the X-ray sources studied bv lFavata et al.l 
((2005!) for the presence of large flares: stars with evidence of very 
long magnetic structures (dashed line), and the complementary 
sample of stars with shorter flaring structures (dotted line). 



compares the abundance distributions of the long-loop 
stars with those characterized by shorter flaring loops, as 
derived bv lFavata et al.l (|2005l ) . While the abundances in 
the short-loop flaring plasma tend to be higher than the 
average, the long-loop objects show systematically lower 
abundance values (also for the other elements not shown 
in figure) . Kolmogorov-Smirnov tests performed between 
these two distributions yield probabilities P — 9% (1%) 
that the Fe (Ne) abundances are drawn from the same 
parent population. These low probabilities suggests that 
we are indeed observing different classes of X-ray sources 
in the two groups, characterized by different chemical 
evolutions or different origins of the flaring plasma. Since 
the statistical significance of this result is not very high, 
specific time- resolved analyses of individual flaring events 
are required to confirm it. 

7. COMPARISON WITH OLDER MAGNETICALLY 
ACTIVE STARS 

We return to the ensemble properties of the X- 
ray bright, count-limited sample of 35 ONC stars 
to compare with other active stars for which high- 
quality X-ray spectroscopy is available. Fig. [12] 
shows the abundances ordered by FIP derived for our 
COUP sources with four comparison stars: the clas- 
sical T Tauri star TWA 5 in the TW Hya associa- 
tion (|Argirofii et alj [20051 ). th e weak-line T Tauri star 
PZ Tel in the (3 Pic asso ciation (lArgiroffi et al.ll2004l). the 
ZAMS star ABDor (jSanz-Forcada. Maggio. fe Miceial 

2003b [), and the active bina ry system V851 Cen 

(|Sanz-Forcada. Favata. fe Miceial (20041) . The abun- 
dances for these four stars were derived from high- 
resolution grating spectra taken with Chandra and/or 
XMM-Newton, and hence with techniques more refined 
than the global spectral fitting approach used here. 

The similarity of the abundance patterns vs. FIP in the 
Orion and older stars is striking, except for the discrepant 
calcium abundances. The X-ray bright ONC stars share 
with other active stars a characteristi c Ne/Fe abundance 
ratio severa l times the solar ratio ([Drake et alJ 120011 : 
iGiidell [20041 ). This behaviour is often attributed to an 
underabundance of low-FIP elements with respect to 
high-FI P elements, i.e. to t he so called "inverse FIP 
effect" ([Brinkman et alJfeOOll ). Such a behaviour was re- 
cently observed also in you ng, weak-line T Tauri stars 
([Argirofli et al.l l2004l . l2005h . and now we find it in the 
X-ray bright ONC stars. 

In recent years, the number of coronal sources 
with available abundance determinations has increased 
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Fig. 12. — Comparison of abundances derived for the COUP 
sources (box plots of the fitting results for the count-limited sub- 
sample) with the abundances obtained from the analysis of high- 
resolution grating spectra of four active stars: TWA5, PZ Tel, 
ABDor, and V851 Cen. 



steadly, but we are far from a clear assessment of the 
phcnomclogy. In fact, a variety of abundance patterns 
has been observed, with more or less pronounced devi- 
ations from both the classical solar FIP effect and the 
inverse FIP effect in its original version. For example, 
the four comparison stars in Fig. rT2] show the inverse- 
FIP effect between iron and neon but they are all char- 
acterized by relatively high abundances for the low-FIP 
elements Ca, Ni, and Mg. Some star-to-star differences 
in the abundance pat terns appear to be l inked to the 
stellar activity level (lAudard et al.l l2003t IGiidell l2004t 
iGarcia- Alvarez et al.l I2006T ) , but again with striking ex- 
ceptions. iWood fe Linskvl (12006) have recently reported 
the case of the binary 70 Oph where the primary shows a 
prominent solar-like FIP effect while the secondary has 
no FIP bias or possibly a weak inverse FIP effect, in 
spite of the similarity between the two stars in all other 
respects. 

8. REVISED TREATMENT OF STANDARD 
ABUNDANCES 

Part of the confusion is likely due to our ignorance of 
true stellar photospheric abundances of the magnetically 
active stars, which are usually assumed to be solar, or at 
least with the same ratios as in the solar photosphere. 
When proper stellar abundance measurements are em- 
ploye d, the abundance vs. FIP pattern is no longer very 
clear ([Sanz-Forcada et al.ll2004f ). 

8.1. Orion photospheric and nebular abundances 

In the case of the ONC, an assessment of the pho- 
tospheric composition is available only for a handful of 
stars. The che mical evolution of t he Orion as sociation 
was studied bv lCunha fe LambertJ (| 19921 11994[ ) who de- 
rived photospheric CNO, Si and Fe abundances for 18 
main-sequence B stars. One of the results of these early 
works was that the spread in O and Si abundances was 
larger than expected based on the measurement uncer- 
tainties. This is thought to indicate a real spread among 
stars of different ages, caused by self-enrichment of the 
nebula as a consequence of supernova explosions within 
the Orion association. However, the result could also be 
affected by systematic errors in the analyses, due to ap- 
proximations in the adopted NLTE model atmospheres 
and line-blanketing effects. In fact, detailed calculations 
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of the oxyge n abundances for 3 B star s in Orion recently 
presented by ISimon-Dfaz et"aH (|2006l) yield values lower 
by ~ .2 dex with respect to those of lCunha fe Lambert! 
(1 19941) for the s ame s tars. 

ICunha et "all (|2006f ) report NLTE Ne abundances for 
11 B-type stellar members of Orion and found a ho- 
mogeneous abundance of neon, A(Ne) = 8.27 ± 0.05, 
and oxygen, A(0) = 8.70 ± 0.09 (in a log scale where 
A(H) — 12). For the same sample, we have computed 
av erage Si and Fe abundan ces from the measurements 
of lCunha fe Lambert! (fl99l ): A(Si) = 7.16 ± 0.15, and 
A(Fe) = 7.46 ± 0.12. These values will be used in the 
next section. 

For later-type stars, ICunha et al.l (|1998l ) determined 
NLTE oxygen and LTE Fe abundances from optical spec- 
troscopy of 9 pre-mai n-sequence F and G Orion stars. 
ICunha fe Smith! (|2005h report a study of fluorine, C and 
O abundances in 3 Orion K-M dwarfs. These works in- 
dicate that the solar-type stars of the Orion association 
all have the same Fe abundance: A(Fe) = 8.40 ± 0.15. 
In contrast, the oxygen abundance appears to vary from 
star to star with a large spread (^4(0) = 8.94 ± 0.36 for 
the full sample) , and we will not consider these measure- 
ments in the following. 

Finally, we consider the composition of the Orion Neb- 
ula, which is the brightest and neares t Galactic H II 
region in the sky. lEsteban et alJ (|2004f ) present echelle 
spectrophotometry of a region S-W of 8 1 Ori C and de- 
rive abundances of several ionic species, including Ne I 
and Ne II, from collisionally-excited lines or recombina- 
tion lines. Their analysis takes into account spatial vari- 
ations of the temperature structure of the nebula and, 
applying ionization and dust-depletion corrections, ob- 
tain abundances of several gas-phase elements. 

8.2. Revised solar abundances 



lAsplund. Grevesse. fe Sauval {2005) present a detailed 
3-D hydrodynamic modeling of the solar atmosphere and 
find that the abundances of many elements (including 
C, N, O, Ar, Ne, and Fe) need to be revised down- 
ward by factors of 1.5-2.4 from the widely-used compi- 
lation of lAnders fe Grevesse! |l989). But the new so- 
lar composition implies lower opacities and produces 
a severe inconsistency between the standard solar in- 
terior model and pre c ise helioseismology measurements 
(jAntia fe Basul l2005t iBahcall et ail |2005|) . One solu- 
tion to this conundrum is to revise upward the poorly- 
known Ne abundance in the Sun, bringing it closer to 
the values often measured in stellar coronae. From 
high-resolu tion X-ray spectroscop y of several active late- 
type stars, iDrake fe Testa! ([2005) showed that the coro- 
nal Ne/O abundance ratio is, on average, a factor 
2.7 times high e r tha n the solar value recommended by 
Asplund et al. (2005|). A similar result was obtained by 
Cunha et alJ (|2006f l. who presented measurements of the 
photospheric Ne abundance in a sample of B-type stars 
in Orion, and obtained a Ne/O abundance ratio a factor 
2.5 higher than the most recent solar value. 

We also find the COUP sources in the count-limited 
sample we have obtained a median Ne/O ratio of 0.33, i.e. 
a facto r 2.2 higher than the lAsphmd. Grevesse. fe Sauvall 
(2005) value. However, this ratio suffers large scatter in 
individual objects (the central 68% of the data span the 
range 0.11-1.18), possibly because oxygen measurements 



are the most affected by uncertainties in the amount of 
absorption and in the amount of low-temperature plasma 
(see Appendix A). 

8.3. ONC coronal abundances with revised standard 
abundances 

The boxes in Fig. [13] shows our Orion coronal abun- 
dances inferred from th e COUP X-ray spec t ra wit h re- 
spect t o the traditional lAnders fe Grevesse! (|1989j ) and 
revised lAsplund et alJ (|2005f ) solar abundances. We have 
not adjusted th e solar neon abundan ce as suggested 
by the work of IDrake fe Testa! (|2005f) . The boxplot 
shows that the inverse-FIP abundance pattern for our 
ONC stars is still pr esent, and is e ven s lightly more 
pronounced, with the lAsplund et al.l (|2005l ) solar abun- 
dances. 

The points with error bars in Figure [13] show average 
stellar photospheric abundances and nebular abundances 
(not corrected for d ust-locking effects) a s described in 
§ 18.11 scaled by the lAsplund et alJ (|2005l) revised solar 
abun dances (a similar scaling is made by (jEsteban et al.l 
12004) ). Our COUP coronal abundances for the high-FIP 
elements S, O, Ar, and Ne are very similar to those of the 
nebula, and also show good agreement with the stellar 
photospheric values for Si, O, and Ne. Thus, while we 
do find a strong inverse FIP effect with respect to solar 
elemental abundances, the effect disappears when Orion 
photospheric and nebular abundances are considered. 

However, discrepancies are found in the iron abun- 
dances which appear significantly lower in the X-ray 
coronal plasma than in the stellar photospheres. The 
very low Fe abundance found for the gaseous nebula can 
be attributed to heavy depletion into grains. If the value 
derived for the B-type and F-G stars is indeed repre- 
sentative of the iron abundance in the photospheres of 
all the late-type ONC stars, Figure [TBI suggests that the 
coronae of these stars are depleted in iron by a factor 
1.5-3. 

9. SUMMARY AND CONCLUSIONS 

• The coronal temperatures and elemental abun- 
dance pattern in X-ray luminous ONC stars is re- 
markably similar to that found from the analysis 
of high-resolution grating spectra of older magnet- 
ically active stars. Hence accretion or the presence 
of circumstellar disks does not appear to affect the 
X-ray production mechanism or plasma. The abun- 
dance of calcium is a possible exception: it appears 
to be extremely low in about 70% of the ONC stars 
we have studied. However, it is difficult to reliably 
measure low calcium abundances with the available 
CCD spectra. 

• Comparison of the observed abundance distribu- 
tions among different stars with simulated distribu- 
tions indicate that all stars may actually have the 
same abundance values, i.e. the abundance spread 
for each element is compatible with the statistical 
uncertainties. Nonetheless, our results also suggest 
possible systematic differences between the abun- 
dance distributions for selected subsamples (e.g. 
those of the sources with short vs. long flaring 
magnetic structures) , which require a specific time- 
resolved spectral analysis to be confirmed. 
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Fig. 13. — Comparison between the abundances derived for the COUP sources in the count-limited subsample — sc aled by two different 
sets of reference so lar values as i ndicated - with th e abundances obtained for 1 1 B-type stars in th e Trapezium (jCunha & Lambertl ll994l ; 
ICunha et aT1l2006h . 8 F-G stars IjCunha et al.lll998T ). and for the Orion nebula (Eateban et al.ll200l) . 



■ The ensemble properties of the COUP X-ray 
brightest sources confirm the low metallicity of 
the coronal plasma with respect to the solar pho- 
tospheric val ue: the median Fe abund ance is w 
0.2 times the lAnders fe Grevessel |l989) value, or 
~ 0.3 times t he m ost recent determination by 
lAsplund etai] (|2005t ). At the same time, the Ne/Fe 
abundance ratio is significantly higher than the so- 
lar one, with a median value 5-7, depending on 
the assumed set of solar abundances. 
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• The X-ray brightest COUP sources show a clear 
pattern of abundances vs. FIP. Extensive simula- 
tions make us confident about the robustness of 
this result. If the solar photospheric abundances 
are adopted for reference, the low FIP elements 
(Mg, Fe, and Si) appear to have similar low abun- 
dances, 0.2-0.3 times the solar values, while Ni 
and the high-FIP elements (S, O, Ar, and Ne) 
appear to have higher abundances, 1-2 times the 
solar ones. However, comparison with abundance 
measurements obtained by means of optical spec- 
troscopy of members of the Orion association indi- 
cates a good agreement between photospheric and 
coronal abundances for Si, O, and Ne, while iron is 
significantly depleted in the X-ray emitting plasma 
with respect to the stellar photospheres, by about a 
factor of 3. We conclude that there is no clear FIP- 
related behaviour of the hot plasma abundances in 
the X-ray bright Orion stars, when proper stellar 
photospheric abundances are taken into account. 
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Fig. 14. — Box plots of the best-fit temperatures (left) and normalizations (right) derived by fitting 1000 simulated 3-T model spectra, 
each with 36,000 (within Poisson statistics). The upper and lower edges of each box comprise the central 68% of the data, the central value 
is the median. Diamonds mark the values in the input model. 



APPENDIX 

SIMULATIONS TO VALIDATE SPECTRAL MODELING 

We have performed several simulations, each with 1000 realizations, based on known 2-T, 3-T, or multi-temperature 
input models and different elemental abundance distributions. They are employed to assess realistic scatter of best-fit 
parameters, in particular on the abundances, and to test for bias in the fitted parameters. We have also performed 
specific simulations in which the photon counting statistics of the observed spectra is taken into account. The selected 
simulations are tailored for comparison with results based on our count-limited sample with more than 10 total 
extracted counts. All simulations use the same ancillary response file (i.e. instrument effective area) and response 
matrix belonging to a real COUP source near the center of the ACIS field of view. The simulated spectra were rebinncd 
as the actual data (§[3]) and the spectral fitting was performed on the same fixed energy range (0.5-8 keV) with the 
same XSPEC procedures. The background spectrum associated with the same source was used as a template in all 
the cases; it contributes 0.3% of the total source+background counts and thus has negligible effect on the results. 
In all simulations, we assume hydrogen column density in the ISM photoelectric absorption model component of 
Ah = 3x 10 21 cm -2 , which is near the median value found for the COUP sources in our low-absorption sample (§ 14. 1| . 

The simulations described below are sorted by increasing complexity, so to explore different sources of uncertainty. 
For each simulation, we state the issue we have tested and we show the relevant results. These simulations validate 
both our ability to recover the correct abundance pattern from the analysis of ACIS spectra, and that the observed 
pattern does not arise in a spurious fashion by our analysis process. 

The first simulation assumes a 3-T input model having temperatures, emission measures, and abundances set to 
values near the median of the distributions obtained for our sources in the low-absorption subsample. The three 
plasma components have kT\ — 0.4 keV, kTi — 0.8 keV, and kT$ = 2.6 keV (logT = 6.7, 7.0, and 7.5 K), and emission 
measure ratios EMz/EMx = 1.5 and EM3/ EM\ = 3.5. All simulated spectra here have 36,000 counts before applying 
Poisson noise, which is near the median value of the spectra which required 3-T best-fit models. Elemental abundances 
were fixed to the median values determined by fitting the sources in our count-limited sample; the Ca abundance, in 
particular, was set to zero. 

Figure [14] shows the distributions of temperatures and volume emission measures derived by fitting the simulated 
spectra, while Fig. [15] shows the distributions of the abundances, sorted by First Ionization Potential (FLP) of the 
relevant elements. The boxplots indicate the range covered by the central 68% of the values. Since the simulation 
is based on a perfect alignment of the input and fitted spectral model, the scatter in parameter values serves as a 
reference for "the best we can do" with Chandra/ ACIS spectra. 

Figure [151 also shows the distributions of abundances derived from simulations of 3-T spectra with 16,000 counts, 
which is the median for all the sources in our count-limited sample. As expected, the width of the distributions is 
slightly larger for all elements than seen in simulations with 36,000 counts. Simulations based on 2-T rather than 3-T 
input spectral models and simulations using sources with the same distribution of counts as in our sample give results 
very similar to those in Figure [151 little additional spread in the derived abundance distributions is introduced by the 
different source model spectrum or by the photon counting statistics of the real COUP dataset. 

In all cases, we see very little bias in the derived spectral parameters; that is, the median values of the distributions 
lie close to (usually within ±10% of) the input values. Temperature estimates become increasingly inaccurate for 
the lower temperature components (logT < 7.0 K, which also have associated emission measures lower than for the 
high-temperature component). Considering the results of all simulations, we find that Fe, Si and S abundances are 
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Fig. 15. — Box plots of the best-fit abundances derived by fitting 1000 simulated 3-T model spectra for sources with 36,000 and 16,000 
counts to illustrate the effects of signal strength. Diamonds mark the values in the input model, and the meaning of the boxes is the same 
as in Fig. ITU 
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Fig. 16. — Box plots of the best-fit abundances derived by fitting 1000 simulated 2-T model spectra with different input abunda nce 
patterns: (left) 0.2x solar abundances; and (right) l.Ox solar abundances with low Fe abundance. Symbols and other details as in Fig. 1141 



generally accurate within 40-80% relative errors, while Ni, O, Ar and Ne have somewhat lower accuracy by factors 
1.8-2.8; Mg shows the largest scatter, with uncertainties up to a factor 10 in simulations with 2-T models and low 
(w 14, 500 counts) photon counting statistics; in the case of the Ca, whose input abundance was assumed to be zero, 
the statistical fluctuations make the best-fit result the most uncertain with any value between and 0.8 acceptable. 

For the cases of Fe and Ne, we have verified that the uncertainties indicated by the simulations are slightly larger 
(by factors 1.2-1.4, on average) than the XSPEC errors, computed at the 90% confidence level for single parameter. 
However, the results presented here show that the uncertainties on the best-fit abundances for all elements are suf- 
ficiently small to recover the input abundance pattern vs. FIP. For example, we ascertain that the Ne abundances 
exceed those of Fe with very high degree of confidence. The scatter on abundance ratios with iron is even lower than 
on individual abundances, because all abundance measurements are correlated with the iron one to a certain degree 4 : 
in fact, the apparent abundance ratio Ne/Fe~6 is affected by « 30% uncertainty, according to our simulations. 

Figure [16] explores the possibility that the interactions between parameters might create the observed FIP abundance 
effect as a form of systemic bias in the fitting procedure. Here we show the abundances emerging from 2-T input models 



4 This is due to the common inverse proportionality between abundances and plasma emission measure while the source count rate is 
fixed. 
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Fig. 17. — Boxplots of temperatures, normalizations, H column densities, and abundances obtained by simulating 1000 spectra with 
an underlying V-shaped emission measure distribution vs. temperature, and fitted with 3-T models. In each panel, diamonds show the 
pattern of input parameter values. 



assuming 0.2 solar abundances for all elements. No trend is evident in the reconstructed abundance pattern. A similar 
result is found for input models with 1.0 solar abundances, exhibiting less scatter due to the stronger emission lines. 
Note, in particular, that an high abundance of Ca could be correctly determined if it were present, but little could be 
inferred for the Ni abundance in this situation (quite different from the case of the COUP sources), due to the strong 
blends with lines from the other elements. 

Our final simulations test our ability to model more complex thermal distributions. Simple 2-T or 3-T models are only 
approximations to the ac tual thermal structurin g of real coronae that must have continuous distributions of emission 
measure vs. temperature. iSchmitt fc Nessl IpOOl warn that coronal abundances may be inaccurately estimated without 
continuous temperature distribution models, but this caveat applies to the analysis of high-resolution grating spectra. 
Fig. [T7| shows a representative simulation in which we assumed a V-shaped emission measure distribution over the 
temperature range logT = 6.0-7.6 K, with a minimum around logT = 6.5 K, resembling distributions determined 
from high-resolution X-ray spectroscopy of very active stars (e.g., Sanz-Forcada, Brickhouse & Dupree 2003). The 
results of this simulation suggest that, although the 3-T model can not effectively recover the actual emission measure 
distribution, the input pattern of abundances can still be reliably measured. Note, in particular, that the inferred Ne 
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is sometimes underestimated but is not significantly overestimated with respect to the input value. Distributions more 
skewed toward low values are those of the O and Mg abundances, but the median is always quite close to the input 
value. 
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8.10 
9.65 
9.46 
9.33 
11.37 
10.84 
10.38 
10.21 
11.20 
9.97 
11.19 
12.45 
10.82 
11.64 
9.65 
11.28 
9.42 
10.61 
10.60 
10.15 
10.52 
9.51 
10.44 
11.19 
10.68 
10.05 
9.66 
10.39 
8.83 
11.30 
8.86 
10.34 
9.91 
9.26 
9.59 
10.61 
10.18 
9.53 
10.06 
9.31 
9.80 
9.17 
11.02 
10.77 
10.06 
9.19 
9.18 
11.67 
8.56 
10.83 
12.27 
9.63 
9.77 
10.68 
11.08 
10.08 
11.07 
11.03 
9.95 
9.64 
11.56 
10.50 
11.12 
11.63 
9.54 
10.17 
9.15 
11.63 
11.76 
10.18 
10.65 
10.67 
11.11 



7.95 
9.46 
8.60 
9.09 
11.05 
10.53 
10.08 
9.53 
10.63 
9.62 
11.01 
12.05 
10.54 
11.29 
9.37 
10.91 
8.99 
10.29 
10.30 
9.86 
10.10 
9.22 
10.12 
10.90 
10.38 
9.87 
9.39 
9.94 
8.26 
11.04 
8.63 
10.08 
9.60 
8.62 
8.78 
10.06 
9.26 
9.14 
9.76 
8.66 
9.43 
8.73 
10.74 
10.32 
9.48 
8.61 
8.89 
11.30 
8.19 
10.60 
11.95 
9.77 
9.82 
10.31 
10.95 
9.72 
10.66 
10.64 
9.68 
9.40 
10.96 
10.19 
10.79 
11.27 
9.41 
9.94 
8.97 
11.37 
11.16 
9.96 
10.39 
10.35 
10.90 



10.34 



9.25 
8.90 



8.87 

8.13 
9.34 
9.66 
9.16 
9.46 
7.59 
9.25 
8.45 



8.62 
7.97 
8.84 
10.59 
7.95 



9.92 

9.24 
10.34 

9.38 



10.93 
9.92 

9.88 
10.60 



Maggio et al. 

TABLE 1 — Continued 



COUP M log* A v A(I-K) EW(Ca) V I J H K s L 

(mag) 



ID 


Sp Type 




(yr) 


(mag) 


(mag) 


(A) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


1433 














17.29 


13.89 


12.02 


11.22 


10.93 


1443 














14.37 


12.48 


11.11 


10.31 


10.08 


1449 














17.63 


14.77 


12.42 


11.35 


10.91 


1463 


K8c-Ml: 


0.59 


5.83 


0.56 


0.50 




14.12 


12.19 


10.78 


10.00 


9.50 


1487 


Ml (SB2) 


0.40 


5.86 


1.65 






15.88 


13.26 


11.55 


10.61 


10.28 


1489 


F9-K0 


2.59 


6.46 


2.06 


0.05 




11.72 


10.30 


9.31 


8.61 


8.40 


1492 


Ml. 5 


0.37 


6.03 


0.11 


0.46 


1.6 


15.13 


13.02 


11.63 


10.88 


10.57 


1499 














15.84 


13.79 


11.39 


10.46 


9.85 


1516 


K1-K4 


1.40 


6.82 


1.31 


-0.20 


1.8 


14.37 


12.77 


11.72 


11.05 


10.89 


1521 


K4 


1.40 


6.67 


1.08 


1.01 




14.31 


12.69 


11.24 


10.36 


9.83 


1568 


KO-Kle 


2.55 


6.34 


0.59 


0.19 




11.30 


10.20 


9.36 


8.84 


8.63 


1595 


M2.5 


0.29 


5.99 


0.00 


0.11 




15.57 


13.25 


11.89 


11.19 


10.97 


1608 


M0.5e 


0.48 


6.23 


0.93 


1.41 


-1.3 


16.01 


13.77 


11.96 


11.06 


10.41 



High-absorption sample 



90 


M0 


0.52 


5.93 


4.97 


0.08 


1.6 


19.09 


15.36 


12.68 


11.44 


10.97 




115 


K7 


0.71 


6.21 


3.83 


0.86 


1.4 


17.99 


14.91 


12.20 


10.91 


10.43 




131 


K5 


1.20 


6.34 


3.95 


0.53 


1.4 


17.12 


14.27 


11.98 


10.95 


10.24 




183 


G: 










0.0 


18.03 


15.68 


12.03 


10.31 


9.23 




223 


K5 


1.19 


6.08 


4.66 


1.04 


1.7 


17.35 


14.22 


11.53 


10.10 


9.34 




262 


K5 


1.13 


6.78 


3.77 


2.24 


2.3 


17.69 


14.91 


11.66 


10.07 


9.30 


8.59 


310 












8.9 


17.91 


14.65 


11.76 


10.40 


9.62 




323 


K6-M0 


0.57 


6.21 


3.86 


1.22 


2.0 


18.50 


15.24 


12.44 


11.11 


10.46 




331 
















16.23 


12.48 


10.55 


9.36 




342 
















14.89 


11.73 


10.25 


9.62 


9.04 


365 


K4-K7 


0.72 


7.79 


0.00 




0.0 




14.37 


11.14 


9.66 


8.82 


7.66 


449 


K7 










0.0 


17.00 


15.38 


12.26 


10.45 


9.38 


8.05 


452 


K0-K3 


2.02 


6.17 


5.36 


0.61 


1.4 


17.04 


13.86 


11.03 


9.73 


8.99 




454 


K2-K7 


2.35 


5.91 


5.85 


0.20 


2.1 


16.85 


13.48 


10.83 


9.61 


9.10 


8.21 


490 


K6-K8 


0.70 


5.62 


4.91 


0.32 


1.2 


17.55 


14.05 


11.40 


10.10 


9.61 




499 


K5-M1 


0.69 


5.92 


2.65 


0.06 


1.2 


16.19 


13.57 


11.68 


10.71 


10.35 




514 














18.20 


15.10 


12.36 


10.93 


10.37 




554 


















12.50 


12.72 


10.42 


8.25 


561 


K5 






0.00 




1.0 




14.58 


11.12 


9.41 


8.34 


6.67 


626 


Ml 


0.41 


6.02 


3.47 


0.62 


0.7 


18.30 


14.97 


12.38 


11.19 


10.78 




649 


M0.5-M2.5 


0.40 


6.03 


4.11 


0.43 


0.0 


19.12 


15.45 


12.64 


11.38 


10.84 


10.46 


655 




















13.40 


10.38 


7.70 


682 


K2-K5 


1.86 


6.05 


2.56 






14.26 


12.12 


10.78 


9.26 


8.63 


7.26 


697 


K5-lateK 










6.1 


15.17 


12.54 


10.25 


9.11 


8.06 


6.84 


707 


M2 


0.33 


5.38 


1.21 


0.83 


1.6 


16.04 


13.34 


11.41 


10.36 


9.77 


8.75 


720 
















16.57 


12.72 


10.93 


10.03 




758 


G5-K0e 


3.00 


6.12 


3.78 


1.50 


-12.3 


13.79 


11.45 


8.78 


7.76 


7.13 


6.16 


766 
















11.00 


9.74 


8.37 


7.35 


5.50 


784 


M1.5-M2e 






0.00 








17.99 


12.20 


11.20 


10.70 


8.63 


874 




















14.31 


12.46 


10.36 


894 


















13.81 


11.82 


10.95 




915 


















13.63 


11.21 


9.98 




939 














18.67 


14.93 


11.49 


9.97 


8.86 




942 


M0 


0.52 


5.90 


5.04 




0.0 


19.06 


15.30 


11.98 


10.43 


9.67 


9.14 


985 


F8-K0 


2.97 


6.17 


2.06 


1.19 




12.23 


10.56 


8.69 


7.75 


7.37 


6.91 


1028 


K2 


1.60 


6.78 


2.72 


1.45 




15.26 


13.20 


10.93 


9.75 


9.11 


8.53 


1035 


















13.62 


11.25 


10.03 




1040 


















13.99 


11.49 


10.16 




1071 


K7-M0 


0.69 


5.92 


1.57 


1.29 


1.6 


15.09 


12.89 


10.38 


9.26 


8.38 


7.38 


1080 


carlyKc-MO 


1.98 


5.25 


7.77 


1.41 


-16.9 


17.82 


13.48 


9.65 


7.72 


6.43 




1114 


K0-K5: 










-1.5 


15.72 


12.46 


9.80 


8.58 


8.04 




1140 














19.59 


15.89 


12.80 


11.28 


10.40 


9.57 


1158 


Ml: 


0.41 


6.08 


2.65 




0.0 


17.67 


14.66 


11.90 


10.54 


9.70 


8.51 


1161 


M0 


0.54 


6.43 


0.00 


1.12 


1.5 


15.21 


13.44 


11.80 


10.90 


10.51 


10.11 


1304 


























1309 


K1-K4 


1.90 


6.44 


4.33 






16.37 


13.64 


10.64 


9.39 


8.75 


7.87 


1335 


K8c 


0.64 


6.64 


1.10 


2.18 


-2.5 


16.32 


14.18 


11.97 


10.74 


9.99 


8.65 


1341 
















14.42 


11.42 


9.96 


9.32 


8.78 


1343 














17.65 


14.63 


12.02 


10.61 


9.78 


8.83 


1354 


G0-M3 


2.34 


6.55 


1.70 


0.18 




11.97 


10.61 


9.67 


9.04 


8.80 


8.25 


1380 


K4 


1.52 


6.27 


3.78 


0.87 




16.28 


13.61 


11.69 


10.31 


9.27 


8.25 


1382 


M0 


0.52 


5.88 


1.97 


1.25 


0.0 


15.94 


13.38 


11.31 


10.17 


9.45 




1391 


Ml 


0.41 


6.03 


3.80 


1.53 


1.4 


18.69 


15.23 


12.14 


10.57 


9.92 




1409 


K6-K8e 


0.74 


7.10 


0.00 


3.01 


-6.3 


15.36 


13.98 


11.70 


10.15 


9.20 


8.07 


1410 


Ml 


0.36 


7.56 


0.57 


2.30 


0.0 


18.54 


16.34 


13.60 


12.31 


11.85 




1421 


M0 


0.54 


6.13 


1.04 


1.10 


0.9 


15.63 


13.43 


11.60 


10.71 


10.33 




1444 


K8c 


0.60 


6.19 


1.10 


0.73 


-4.1 


15.55 


13.41 


11.78 


10.94 


10.43 




1456 
















16.83 


13.36 


11.44 


10.57 




1462 


Kl 


2.54 


6.25 


3.09 


0.18 




13.95 


11.82 


10.12 


9.25 


8.92 





Coronal abundances in Orion 
TABLE 1 — Continued 



COUP 
ID 


Sp Type 


M 
(M e ) 


log* 

(y) 


A v 
(mag) 


A(I-K) 

(mag) 


EW(Ca) 
(A) 


V 
(mag) 


/ J 

(mag) (mag) 


H K s 
(mag) (mag) 


L 

(mag) 


1466 


K3-K5 


1.40 


6.65 


1.96 


0.42 


1.4 


15.15 


13.19 11.69 


10.78 10.38 





TABLE 2 

Spectral analysis results 



COUP AT H 



T 2 



T 3 



ID 




keV 


keV 


keV 


7 


1.5 


0.2 


0.7 


2.2 


9 


2.6 




0.8 


3.1 


11 


3.9 




0.6 


5.4 


23 


2.4 


0.1 


0.8 


2.4 


27 


1.7 




0.7 


2.7 


28 


1.5 




0.6 


3.4 


43 


5.3 




0.3 


2.6 


62 


2.6 




0.6 


3.1 


66 


1.7 




0.7 


3.2 


67 


2.0 




0.7 


2.4 


71 


0.0 




0.7 


2.9 


101 


1.6 




0.9 


3.7 


108 


3.1 




0.4 


2.9 


112 


3.2 




0.6 


2.5 


113 


2.9 




0.8 


2.3 


139 


1.6 




0.8 


2.7 


141 


1.8 




0.8 


4.0 


150 


2.8 




0.7 


2.8 


152 


0.5 




0.7 


2.7 


173 


1.9 




0.8 


2.7 


177 


4.6 




0.8 


3.1 


188 


2.7 




0.7 


2.5 


202 


2.4 




0.5 


2.2 


205 


1.9 




0.8 


4.0 


270 


1.3 




0.8 


2.3 


328 


1.5 




0.6 


2.0 


343 


1.6 




0.7 


3.1 


382 


2.7 




0.6 


2.3 


387 


2.4 




0.7 


2.5 


417 


1.9 




0.8 


3.0 


431 


5.1 




0.9 


2.4 


459 


0.2 




0.7 


2.7 


470 


1.4 




0.8 


2.7 


567 


2.1 




0.7 


3.0 


579 


4.5 




0.4 


4.4 


597 


0.0 




0.7 


2.4 


600 


2.7 




0.6 


2.4 


648 


3.8 




0.6 


2.7 


669 


3.6 




0.6 


3.0 


670 


4.4 


0.5 


2.2 


7.7 


672 


2.6 




0.4 


2.5 


718 


3.8 




0.6 


2.9 


752 


2.8 


0.4 


1.3 


4.2 


753 


5.3 




0.4 


2.1 


761 


5.5 




0.8 


2.8 


801 


5.9 




0.4 


3.6 


828 


3.5 




0.5 


4.0 


848 


2.7 


0.5 


1.7 


10.8 


867 


4.1 




0.6 


2.7 


945 


1.0 




0.4 


1.9 


960 


1.1 




0.7 


3.3 


971 


0.3 


0.2 


0.8 


2.6 


982 


2.7 




1.3 


4.7 


997 


3.6 




0.7 


3.0 


1002 


0.6 




0.4 


1.9 


1083 


5.0 




0.6 


2.6 


1111 


3.8 




0.6 


2.8 


1127 


5.4 




0.5 


2.0 


1143 


1.8 




0.4 


2.2 


1151 


1.6 




0.7 


2.6 


1246 


4.3 




0.4 


3.4 


1248 


3.0 




0.8 


2.5 


1252 


4.7 




0.7 


3.3 


1261 


2.8 




0.6 


3.5 



fcl &2 
(b) (b) 



k 3 
(b) 



O Ne 



Abundances (solar units) 
Mg Si S Ar Ca 



Fe Ni 



DoF (c) 



0.6 



2.2 
1.6 
0.9 
3.6 
0.3 
0.5 
18.6 
1.0 
0.2 
0.4 
0.6 
0.4 
0.2 
0.2 
1.1 
0.6 
1.0 
0.3 
0.3 
1.0 
1.2 
1.8 
0.3 
1.1 
0.3 
1.2 
1.1 
0.2 
0.9 
0.5 
1.0 
0.1 
0.5 
0.5 
0.1 
0.2 
0.2 
0.7 
0.5 
2.9 
0.5 
0.7 
1.2 
0.7 
0.3 
0.6 
0.3 
0.6 
0.4 
0.3 
0.0 
1.2 
1.2 
0.7 
0.9 
0.7 
0.2 
0.4 
0.4 
1.5 
0.3 
0.8 
0.3 
0.1 



3.8 
2.8 
0.7 
5.6 
0.5 
2.5 
1.0 
0.9 
0.6 
0.7 
0.4 
1.3 
1.0 
0.7 
2.0 
0.4 
1.0 
0.6 
0.4 
0.8 
0.4 
3.4 
0.6 
0.5 
0.4 
1.0 
2.8 
0.3 
2.1 
0.6 
2.7 
0.6 
0.9 
1.0 
1.0 
0.4 
0.4 
2.8 
2.4 
1.0 
1.2 
2.5 
3.8 
1.0 
0.8 
3.6 
4.2 
0.3 
1.9 
0.5 
0.4 
3.6 
1.2 
2.0 
0.7 
2.1 
0.8 
1.1 
1.7 
2.1 
1.0 
1.8 
0.6 
0.5 



0.4 
0.9 
0.1 
0.8 
0.7 
0.6 
0.0 
0.4 
1.2 
0.9 
0.3 
0.9 
0.7 
1.2 
1.1 
0.5 
0.8 
0.9 
0.7 
1.0 
0.0 
0.6 
0.3 
0.3 
1.4 
0.5 
1.2 
1.9 
0.5 
0.5 
0.9 
0.7 
1.5 
0.7 
0.1 
1.0 
0.8 
1.0 
1.0 
0.6 
0.3 
0.0 
0.3 
0.5 
7.3 
4.8 
0.3 
0.5 
0.3 
0.4 
1.4 
0.6 
0.0 
0.0 
0.2 
0.3 
2.6 
0.0 
0.3 
0.5 
1.3 
0.0 
1.4 
2.2 



1.0 
0.6 
0.3 
0.9 
1.4 
1.1 
0.0 
0.2 
1.3 
1.1 
0.1 
0.3 
2.1 
2.1 
0.9 
0.2 
0.5 
1.1 
0.4 
0.1 
0.0 
0.9 
0.7 
0.2 
0.8 
0.8 
1.7 
2.3 
1.1 
0.7 
1.4 
2.4 
0.6 
1.3 
2.5 
1.6 
1.6 
2.3 
2.3 
2.0 
1.4 
1.2 
0.9 
0.9 
2.1 
4.6 
4.4 
0.7 
2.3 
1.4 
2.3 
1.6 
1.4 
1.0 
0.8 
1.5 
2.4 
1.3 
1.8 
0.8 
1.8 
0.7 
1.5 
1.9 



0.3 
0.4 
0.2 
0.2 
0.2 
0.2 
0.0 
0.1 
0.0 
0.5 
0.0 
0.3 
0.9 
0.4 
0.1 
0.0 
0.1 
0.2 
0.0 
0.0 
0.1 
0.1 
0.0 
0.0 
0.2 
0.2 
0.1 
0.7 
0.2 
0.3 
0.3 
0.8 
0.2 
0.5 
0.0 
0.0 
0.8 
0.4 
0.5 
0.6 
0.7 
0.2 
0.2 
0.5 
0.5 
1.5 
0.4 
0.0 
0.6 
0.7 
0.4 
0.0 
0.0 
0.3 
0.1 
0.1 
0.7 
0.4 
0.5 
0.0 
0.2 
0.1 
0.6 
0.5 



0.2 
0.4 
0.7 
0.1 
0.2 
0.3 
0.1 
0.2 
0.0 
0.2 
0.0 
0.0 
0.5 
0.5 
0.1 
0.0 
0.2 
0.2 
0.1 
0.1 
0.1 
0.1 
0.2 
0.1 
0.4 
0.1 
0.3 
1.0 
0.2 
0.1 
0.4 
0.7 
0.2 
0.4 
1.1 
0.5 
0.9 
0.5 
0.5 
0.6 
0.8 
0.3 
0.1 
0.3 
0.3 
1.3 
1.2 
0.3 
0.7 
0.4 
0.4 
0.4 
0.5 
0.3 
0.2 
0.2 
0.5 
0.5 
0.4 
0.1 
0.5 
0.2 
0.6 
0.6 



0.3 
0.6 
2.9 
0.5 
0.6 
0.9 
0.6 
1.0 
1.5 
0.3 
0.9 
0.8 
0.0 
0.6 
0.8 
0.4 
1.0 
0.9 
0.6 
0.7 
0.3 
0.4 
0.4 
1.0 
0.6 
0.4 
1.2 
1.7 
0.7 
1.1 
0.6 
0.4 
1.1 
0.8 
1.6 
1.1 
1.7 
1.2 
0.8 
1.0 
1.3 
0.8 
0.3 
0.5 
0.8 
2.0 
1.6 
0.9 
0.6 
0.9 
1.2 
0.5 
0.7 
1.3 
0.5 
0.7 
1.0 
0.5 
0.6 
0.9 
0.1 
0.7 
1.1 
1.2 



0.7 
1.3 
2.5 
0.6 
1.3 
0.5 
0.0 
1.6 
0.6 
1.4 
0.9 
1.4 
0.6 
0.9 
0.6 
1.5 
1.1 
1.7 
0.0 
0.9 
1.1 
0.5 
0.9 
2.7 
0.8 
0.5 
0.8 
0.2 
0.4 
1.2 
0.6 
1.6 
1.1 
1.4 
3.8 
1.7 
0.3 
0.9 
0.9 
1.7 
0.5 
0.0 
0.8 
1.0 
1.0 
2.9 
1.4 
1.3 
0.5 
0.0 
2.5 
0.9 
0.4 
0.2 
0.1 
0.8 
0.7 
1.5 
1.0 
0.9 
0.3 
0.7 
0.6 
0.0 



0.2 
0.0 
0.0 
0.1 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
1.8 
1.4 
0.1 
0.9 
0.0 
0.0 
0.0 
0.3 
1.7 
0.3 
0.1 
0.0 
0.0 
0.8 
0.5 
0.0 
0.4 
4.9 
0.3 
1.4 
0.8 
0.0 
0.0 
0.5 
0.0 
5.7 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.6 
0.0 
0.0 
0.0 
0.0 
3.0 
0.0 
0.0 
3.2 
0.0 
0.3 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.1 
0.0 



0.2 
0.3 
0.1 
0.2 
0.2 
0.2 
0.0 
0.0 
0.1 
0.3 
0.1 
0.1 
0.3 
0.6 
0.2 
0.1 
0.2 
0.2 
0.1 
0.2 
0.1 
0.1 
0.1 
0.1 
0.3 
0.1 
0.2 
0.4 
0.2 
0.1 
0.3 
0.4 
0.4 
0.2 
0.2 
0.2 
0.3 
0.2 
0.3 
0.2 
0.2 
0.2 
0.2 
0.2 
0.9 
0.5 
0.4 
0.1 
0.4 
0.2 
1.0 
0.2 
0.0 
0.1 
0.1 
0.2 
0.4 
0.1 
0.1 
0.1 
0.3 
0.1 
0.5 
0.4 



0.7 
1.7 
0.7 
0.7 
0.4 
0.7 
1.6 
1.8 
0.6 
0.8 
1.3 
2.8 
2.9 
1.6 
1.3 
0.0 
0.1 
1.1 
1.5 
0.2 
0.1 
1.1 
2.1 
0.0 
0.0 
1.7 
0.1 
3.9 
1.0 
1.2 
2.8 
1.9 
2.2 
1.5 
0.0 
3.3 
3.1 
1.7 
0.6 
1.2 
1.0 
1.8 
1.5 
0.0 
1.2 
2.4 
1.4 
1.2 
3.7 
1.3 
3.2 
1.9 
0.3 
1.5 
1.3 
1.8 
2.1 
0.1 
1.2 
1.3 
1.3 
1.7 
1.7 
2.7 



1.2 
1.3 
1.6 
1.2 
0.9 
1.2 
1.1 
1.2 
1.2 
1.0 
1.0 
0.8 
0.9 
1.1 
1.3 
0.9 
1.1 
1.3 
1.2 
1.0 
0.8 
1.3 



1.6 
0.9 
1.2 
0.8 
1.1 
0.9 
1.3 
1.0 
1.6 
1.0 
0.9 
1.8 
1.2 
1.6 
1.5 
1.2 
0.7 
1.0 
1.1 
1.3 
1.1 
1.2 
1.2 
0.9 
1.2 
1.3 
1.1 
1.0 
0.9 
1.2 



1.4 



162 

125 
61 

203 
61 

154 
74 
97 
64 
75 
62 

103 
69 
77 
71 
62 

117 
64 
58 
99 
48 

158 
49 
66 
63 
97 

177 
50 

135 
62 

141 
78 
99 
99 
79 
78 
52 

159 

149 

185 
80 

143 

199 
63 

108 

127 

123 
77 

128 
63 
56 

142 
63 

130 
76 

124 
81 
60 

119 

146 
85 

117 
59 
57 



1.5 
1.7 
0.6 
2.6 
0.3 
1.6 
0.4 
0.5 
0.4 
0.3 
0.2 
0.9 
0.6 
0.4 
0.9 
0.2 
0.8 
0.3 
0.2 
0.5 
0.3 
1.5 
0.2 
0.4 
0.2 
0.4 
1.8 
0.2 
1.0 
0.4 
1.3 
0.3 
0.5 
0.6 
0.7 
0.3 
0.2 
1.6 
1.4 
2.3 
0.6 
1.2 
2.8 
0.4 
0.7 
3.4 
3.2 
0.5 
1.0 
0.2 
0.3 
1.8 
1.1 
1.1 
0.2 
1.0 
0.5 
0.4 
0.7 
1.0 
0.6 
0.8 
0.4 
0.4 
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TABLE 2 — Continued 



COUP 


N n 


Ti 


T 2 


T 3 




fc 2 








Abundances (solar units) 










/x 


ID 


fa) 


keV 


keV 


keV 


(b) 


(b) 


(b) 





Nc 


Mg 


Si 


S 


At 


Ca 


Fe 


Ni 


X? 

VI 


DoF 


(c) 


1269 


1.5 


0.4 


0.8 


1.8 


1.1 


1.4 


2.7 


0.4 


1.1 


0.4 


0.1 


0.3 


0.6 


0.0 


0.2 


0.3 


1.1 


136 


0.9 


1311 


2.1 




0.5 


2.2 




0.2 


0.1 


0.5 


2.0 


0.4 


0.7 


1 . 4 


U. i 


1.6 


0.3 


3.2 


1.0 


51 


0.2 


1350 


2.2 


0.5 


0.9 


2.3 


1.0 


1.2 


2.8 


0.6 


0.8 


0.3 


0.3 


0.6 


0.5 


0.1 


0.2 


1.2 


1,1 


145 


1.3 


1355 


1.5 




0.4 


2.0 




0.5 


0.6 


0.1 


0.8 


0.1 


0.3 


0.4 


0.0 


0.0 


0.1 


1.1 


1.2 


58 


0.2 


1374 


4.2 




0.1 


2.6 




0.1 


0.9 


0.0 


2.1 


0.3 


0.2 


0.8 


U. i 


1.0 


0.2 


1.3 


1.3 


56 


0,1 


1384 


1.0 




0.4 


3.2 




0.2 


2.4 


0.7 


2.7 


0.1 


0.6 


0.4 


0.5 


0.8 


0.2 


2.3 


1.1 


160 


1.5 


1412 


1.6 




0.8 


2.2 




0.3 


0.6 


0.7 


0.7 


0.1 


0.3 


1 . 1 


0.3 


1.8 


0.3 


1.2 


1.0 


74 


0.3 


1424 


3.1 




0.5 


2.1 




0.1 


0.9 


0.3 


1.5 


0.2 


0.2 


0.4 


0.2 


0.0 


0.1 


2.2 


1.1 


74 


0.3 


1429 


1.1 




0.4 


2.0 




0.2 


0.6 


0.3 


2.0 


0.0 


0.3 


0.4 


0.1 


0.0 


0.2 


0.6 


1.1 


52 


0.2 


1433 


4.2 




0.4 


2.5 




0.2 


0.6 


3.7 


4.1 


0.0 


0.3 


0.4 


0.2 


0.0 


0.5 


5.0 


0.9 


71 


0,1 


1443 


2.2 




0.7 


3.1 




0.2 


0.8 


1.7 


2.0 


0.1 


0.8 


1.3 


0.0 


0.4 


0,5 


2.3 


1.2 


97 


0.6 


1449 


4.0 




0.4 


3.3 




0.0 


0.8 


0.0 


4.4 


0.3 


1.0 


1.5 


0.8 


0.0 


0.1 


1.7 


1.2 


61 


0.5 


1463 


0.7 




0.4 


2.5 




0.2 


0.8 


0,1 


1.7 


0.2 


0.4 


0,1 


0.2 


0.0 


0.3 


1.9 


1.1 


79 


0.3 


1487 


1.2 




0.4 


2.2 




0.3 


0.6 


1.0 


2.2 


0.8 


0.7 


0.8 


0.8 


0.7 


0.4 


4.6 


1.2 


60 


0.3 


1489 


1.9 


0.1 


0.8 


2.0 


0.4 


0.5 


1.5 


0.1 


1.2 


0.5 


0.1 


0.7 


0.3 


0.3 


0.2 


1.2 


1.1 


107 


0,5 


1492 


1.3 




0.7 


2.4 




0.2 


0.5 


1.2 


2.1 


0.7 


0.6 


0.9 


0.1 


0.0 


0.3 


2.4 


1.0 


67 


0.2 


1499 


5.9 




0.6 


4.0 




0.2 


0.7 


4.0 


2.1 


0.8 


0.9 


2.1 


1.6 


1,1 


0.6 


0.0 


1.0 


65 


0.7 


1516 


0.3 


0.2 


0.6 


1.9 


0.2 


0.2 


0.6 


0.3 


1.5 


0.3 


0.4 


0.2 


1.7 


0.0 


0.2 


1.0 


1.1 


68 


0.2 


1521 


1.5 




0.5 


2.6 




0.4 


1.0 


0.1 


1.7 


0.5 


0.4 


0.7 


0.0 


0.0 


0.1 


0.2 


1.6 


93 


0.5 


1568 


0.6 


0.8 


1.5 


8.1 


1.9 


4.1 


4.2 


0.8 


0.6 


0.2 


0.4 


0.5 


0.7 


0.7 


0.3 


1.2 


2.7 


265 


5,1 


1595 


0.3 




0.7 


2.4 




0.5 


0.3 


0.7 


0.8 


0.2 


0.3 


1.0 


0.0 


4.2 


0.1 


1.6 


1.8 


60 


0.2 


1608 


2.9 


0.2 


1.3 


5.2 


1.2 


1.3 


0.6 


0.1 


0.6 


0.0 


0.2 


0.1 


1.3 


0.0 


0.2 


0.0 


1.2 


92 


0.6 



a Interstellar absorption H column density, in units of 10 21 cm 2 . Volume emission measure for each thermal component, divided by 4ttD 2 
(with D the source distance), in units of 10 — 10 cm -5 . C Absorbed source X-ray flux, in the 2-8 keV (hard) band, in units of 10~ 13 erg cm -2 
s" 1 . 

TABLE 3 

Abundances in ONC X-ray bright stars 



Low-absorption sample Count-limited sample 
FIP median 68% range median 68% range 
Element (eV) (Ajiders_^ i Grcvcssc (198,9) solar abundance units) 



Ca 


6.11 


0.00 


0.00 


0.90 


0.00 


0.00 


0.54 


Ni 


7.64 


1.32 


0.31 


2.51 


1.37 


0.65 


2.20 


Mg 


7.65 


0.25 


0.03 


0.58 


0.26 


0.09 


0.49 


Fe 


7.90 


0.22 


0.12 


0.37 


0.22 


0.14 


0.33 


Si 


8.15 


0.33 


0.14 


0.70 


0.30 


0.13 


0.55 


S 


10.36 


0.75 


0.43 


1.23 


0.74 


0.45 


1.13 





13.62 


0.57 


0.25 


1.26 


0.60 


0.29 


1.00 


Ar 


15.76 


0.76 


0.20 


1.52 


0.76 


0.49 


1.33 


Ne 


21.56 


1.34 


0.59 


2.23 


1.14 


0.65 


2.27 



